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ABSTRACT 



Genus Epinephelus (Perciformes, Epinephelidae), commonly known as groupers, 
are usually difficult in species identification for the lack and/or change of mor- 
phological specialization. In this study, molecular cytogenetic analyses were firstly 
performed to identify the closely related species Epinephelus bruneus and E. moara in 
this genus. The species-specific differences of both fish species showed in karyotype, 
chromosomal distribution of nucleolar organizer regions (NORs) and localization of 
18S rDNA. The heterochromatin (interstitial C-bands) and distribution pattern of 
telomere (TTAGGG) n in E. bruneus revealed the chromosomal rearrangements and 
different karyotypic evolutionary characteristics compared to those in E. moara. The 
cytogenetic data suggested that the lineages of E. bruneus and E. moara were recently 
derived within the genus Epinephelus, and E. moara exhibited more plesiomorphic 
features than E. bruneus. All results confirmed that E. moara, which has long been 
considered a synonym of E. bruneus, is a distinct species in the family Epinephelidae. 
In addition, molecular cytogenetic analyses are useful in species differentiation and 
phylogenetic reconstruction in groupers. 



Subjects Aquaculture, Fisheries and Fish Science, Evolutionary Studies, Genetics, Marine 
Biology, Taxonomy 
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INTRODUCTION 

The family Epinephelidae comprises approximately 163 grouper species in 16 genera 
{Craig, Sadovy de Mitcheson & Heemstra, 2011). These species are of considerable economic 
value, especially in the coastal fisheries of tropical and subtropical areas (Heemstra 
& Randall, 1993). Taxonomic confusion in the Epinephelidae often occurs due to 
similarities of color patterns and ontogenetic changes in color (Heemstra & Randall, 1993; 
Craig, Sadovy de Mitcheson & Heemstra, 2011). Epinephelus bruneus (Bloch 1793) and 
E. moara (Temminck and Schiegel 1842) are two important aquaculture and commercial 
fish species. However, E. moara has long been considered a synonym of E. bruneus 
due to their similarities in coloration and overlapping in geographical distributions 
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(Heemstra & Randall, 1993; Craig, Sadovy de Mitcheson &Heemstra, 2011). Designation 
of correctly identified new species is important not only to the study of phylogenetic 
relationships, but also to the identification of fry and parent in grouper aquaculture. Based 
on morphological characteristics and molecular comparisons, E. moara has been suggested 
as a valid species (Guo et al., 2008; Guo et al, 2009; Liu et al, 2013). The interspecific 
differences between them were identified mainly based on the skeleton system as well as 
the meristic and morphometric characteristics (Guo et al, 2008). Gene differentiation 
(Guo et al, 2009) and mitogenome analyses (Liu et al, 2013) partially provided molecular 
information confirming their taxonomic status. Chromosomes are the carriers of genetic 
information, and chromosomal recombination plays a vital role in genetic diversity 
Therefore, more other evidences are needed to support the hypothesis that E. moara is 
a valid species, such as molecular cytogenetic analyses. 

Chromosomes are hereditary elements of the complete nuclear genome. Molecular 
cytogenetic studies on chromosomes constitute important approaches for characterizing 
species and reconstructing phylogenetic relationships (Galetti, Aguilar & Molina, 2000; 
Ocalewicz, Woznicki & Jankun, 2008; Cioffi, Martins &Bertollo, 2010; Ruiz-Herrera, Farre 
& Robinson, 2012). Karyological features indicate the evolutionary distance between 
species of different taxonomic categories (Dobigny et al, 2004). The nucleolar organizer 
regions (NORs) were particularly significant in chromosomal evolutionary analyses 
(Miller et al, 1976; Fujiwara et al, 1998). Heterochromatin corresponding to C-bandings 
is normally associated with rearrangements, quantitative variation, and formation of 
new karyotypes (Miklos & Gill, 1982; Rocco et al, 2002). Many taxonomic studies were 
based on the variations and polymorphism of the chromosomes containing major 18S 
rDNA (both active and non-active) (Cioffi, Martins & Bertollo, 2010; Britton-Davidian, 
Cazaux & Catalan, 2012), and minor 5S rDNA (Fujiwara et al, 1998; Mazzei et al, 2004) 
by fluorescent in situ hybridization (FISH). Location of telomeric sequence (TTAGGG)„ 
provides direct evidence for cytotaxonomic studies and chromosomal evolution in fishes 
(Sola et al, 2003; Scacchetti et al, 201 1). Therefore, molecular cytogenetic information has 
provided important contributions to the characterization of biodiversity and the evolution 
of ichthyofauna (Jesus et al, 2003; Vicari et al, 2008). The analyses of different methods 
above present a complete karyotypic picture for organisms. 

The cytogenetic information provided by a variety of approaches will allow us to 
more fully explain the taxonomic and evolutionary statuses, and reveal the inherent 
differences of E. bruneus and E. moara. In this study, the karyotypic techniques, including 
Giemsa-staining, Ag-staining, C-banding and localization of 18S rDNA and telomere 
(TTAGGG) n by FISH, were used to investigate the molecular cytogenetic characteristics of 
E. bruneus and E. moara. Cytogentic data provided a better definition of the specific epithet 
for these cryptic species. Molecular cytogenetic analyses were found to be applicable 
in differentiating between closely related species and reconstructing phylogenetic 
relationships in groupers. 
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Figure 1 Specimen of adult E. bruneus and E. moara. Specimen of adult E. bruneus (Bloch, 1973) 
[standard length (LS), 550 mm] and E. moara (Temminck & Schiegel, 1842) (LS, 650 mm) (Xiamen, 
Fujian, China, ML Guo). 



MATERIALS AND METHODS 
Fish collection and identification 

Twenty-five individuals of E. bruneus (standard length, L$, 140-550 mm) and 24 
individuals of E. moara (L$, 188-650 mm) were collected alive from the coastal waters 
of Fujian, China, and reared in laboratory for one week before analyses. Specimen 
identification was based on external coloration (Fig. 1), skeleton and morphological 
characteristics established in previous studies (Guo et ah, 2008). 

Chromosome preparation, karyptyping and staining analyses 

Fishes were injected with colchicine (3 |.ig/g weight, Sigma) for 30 min. Mitotic 
chromosomes were obtained from cell suspensions of anterior kidney after the fishes 
were anesthetized with tricaine methanesulfonate (MS222, 100 mg/L, Sigma), using the 
conventional air-drying method (Ojima, Hitotsumachi &Makino, 1966). Chromosomes 
were stained using Giemsa and classified as metacentric (m), submetacentric (sm), 
subtelocentric (st), or acrocentric (a) based on the arm ratios (Levari, Fredga & Sandberg, 
1964). The nucleolar organizer regions (NORs) were visualized by Ag-staining (Howell 
& Black, 1980). Heterochromatin was identified by C-banding using barium hydroxide 
method (Sumner, 1972). After the acquisition of anterior kidney, tissue samples (mostly 
dorsal muscle) were collected and stored at —80 °C. All experiments were conducted 
in accordance with the guidelines approved by the Institutional Animal Care and Use 
Committee at Xiamen University. 

Chromosomal probes preparation 

Genomic DNA of all specimens was extracted from muscle tissue using the phenol- 
chloroform method (Sambrook, Fritsch & Maniatis, 1989). 18S rDNA and 5S rDNA 
probes for chromosome hybridization were prepared as follows: The partial coding 
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region of 18S rDNA were amplified using the primers 18S rDNA-F (5'-GTAGTCA 
TATGCTTGTCTC-3') and 18S rDNA-R (S'-TCCGCAGGTTCACCTACGGA^O as 
described by White et al. (1990). The coding region of 5S rDNA were obtained using 
the primers 5S rDNA-F (S'-TACGCCCGATCTCGTCCGATC-SO and 5S rDNA-R 
(S'-CAGGCTGGTATGGCCGTAAGC-a') indicated by Martins & Galetti (1999). PCR 
reactions were performed as following: 94 °C for 4 min, followed by 35 cycles of 94 °C 
for 30 s, 54 °C (for 18S rDNA) or 62 °C (for 5S rDNA) for 1 min, and 72 °C for 1 min, 
and a final extension at 72 °C for 5 min. The nucleotide sequences of 18S rDNA and 
5S rDNA were obtained after cloning into the pMD-18T vector (Takara, Japan), and 
subjected to Blastn in NCBI database (http://www.ncbi.nlm.nih.gov). Telomere probes 
for chromosome hybridization were prepared as follows: telomeric repeat sequences 
(TTAGGG)„ were amplified by PCR using (5'-TTAGGG-3')5 and (5'-CCCTAA-3')5 as 
primers (Ijdo etal, 1991). All probes were labeled with biotin-16-dUTP (Roche, Germany) 
by nick translation according to the manufacturer's instructions. 

Fluorescence in situ hybridization (FISH) 

FISH and probe detection were conducted using methods as described previously ( Wang 
et al., 2010). Briefly, avidin-fluorescein isothiocyanate (FITC) (Sigma, USA) was used 
for signal detection of probes 18S rDNA, 5S rDNA and telomere (TTAGGG)„ based 
on the manufacturer's instruction. Chromosomes were counterstained with 1 |Ag/ml A', 
6 / -diamidino-2-phenylin-dole (DAPI) (Roche, USA) in anti-fade solution of 70% glycerol, 
2.5% DABCO [1,4-Diazabicyclo (2.2.2) octan], and lx standard saline concentration 
(SSC) at pH 8.0. Hybridization signals were observed and analyzed under a fluorescence 
microscope Leica DM-400CCD. 

RESULTS 

Karyotypes and banding patterns 

A total of 1 72 metaphases of E. bruneus and 156 metaphases of E. moara were analyzed to 
determine the karyotype structure. All specimens of E. bruneus and E. moara invariably 
showed the same diploid number of chromosomes, 2n = 48. The karyotype formulas 
of E. bruneus and E. moara were 2m + 4sm + 42a, giving a fundamental number (NF) 
equaled to 54 (Figs. 2A and 2B), and 4sm + 44a, NF = 52 (Figs. 2C and 2D), respectively. 
Chromosomes pairs were numbered based on the relative length. The smallest chromo- 
somes pairs No. 24 were submetacentric chromosomes (sm-3 for E. bruneus and sm-2 for 
E. moara, Figs. 2B and 2D). Chromosome pairs No. 9 in length were sm-2 for E. bruneus 
and sm-1 for E. moara. Chromosome pairs No. 2 were metacentric chromosomes m-1 for 
E. bruneus. Other chromosomes were acrocentric (a) chromosomes for both E. bruneus 
and£. moara. 

Active NORs were identified on the terminal position of short arms or sub-centromere 
regions of those biarmed chromosomes. In E. bruneus, five actively transcribed NORs were 
located on the metacentric and submetacentric chromosomes (Fig. 3A). In E. moara, four 
Ag-NORs were found on the submetacentric chromosomes (Fig. 3B). 



Guo et al. (2014), PeerJ, DOI 10.7717/peerj.412 



4/15 



PeerJ 



B 



m 



sm 



a 



18 

i 

fl > 

2 

n 

4 

10 

0/J 

16 

00 

22 



3 

n 

5 

u 

11 

17 

80 

23 



ft* 

6 

12 

H 

18 

It 

24 



7 

Oft 

13 

H 

19 



Ad 

8 

M 
14 

/II 

20 



8 (I 

9 

90 

15 

06 

21 



4 



t 4 4> * 



sm 



10 

1 

• A 

3 



2 

*A 

4 



A«^ RO 

9 10 

A Q rv ft 

15 IS 

M 

21 22 



»A 0 «o id or 

5 6 7 8 

I^n /v-* or>> a a 

11 12 13 14 

OA 0 0 Aft ftf 

17 18 19 20 

tft « • 

23 24 



Figure 2 Chromosome metaphase and corresponding karyotype of E. bruneus and E. moara. Chro- 
mosome metaphase (left, Giemsa staining) and corresponding karyotype (right) of E. bruneus (A and B) 
and E. moara (C and D). Scale bar = 5 |tm. a, acrocentric; m, metacentric; sm, submetacentric. 



The constitutive heterochromatin was observed in the centromeric and/or peri- 
centromeric region of most chromosomes for both E. bruneus and E. moara. And the 
biarmed chromosome pairs with positive Ag-NORs were coinciding with the positive 
heterochromatin C-bandings. While three pairs of acrocentric chromosomes were almost 
indiscernible in both fish species (Figs. 4A and 4B). However, the significant differences 
of heterochromatin were the heterochromatic blocks found in the interstitial region of the 
long arms of one pair of medium-sized acrocentric chromosome in E. bruneus (Fig. 4A). 

Sequences analyses 

Sequences of 18S rDNA (GenBank accession nos. FJ176793 and FJ176794) and 5S rDNA 
(GenBank accession nos. FJ1 76796 and FJ176795) were amplified from genomic DNA of 
E. bruneus and E. moara. Sequence of 18S rDNA contained partial DNA of gene 18S rRNA. 
Partial DNA sequence of 5S rDNA included the encoding and non-transcribed spacer 
(NTS) region for both fish species. The determined sequences were highly conserved. 
The nucleotide similarities of partial 18S rDNA and 5S rDNA were 100% and 99.99%, 
respectively, for both fish species. The phylogenic neighbor-joining (NJ) trees based on 
partial sequences of 18S rDNA and 5S rDNA strongly support the closed relationship of 
E. moara and E. bruneus (high bootstrap values of 92 and 100). And genera of the order 
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Figure 3 Ag-NORs characteristics of E . bruneus and E. moara with silver staining. Ag-NORs char- 
acteristics of E. bruneus (A) and E. moara (B) with silver staining. Thick black arrows indicate the 
chromosomes No. 9 in length, thin black arrows represent the chromosomes No. 24 in length, and hollow 
arrows show the chromosomes No. 1 in length. Scale bar = 5 |im. 




Figure 4 C-banding patterns of E. bruneus and£. moara. C-banding patterns of E. bruneus (A) and 
E. moara (B). Heterochromatic blocks were observed in the interstitial region of the long arms of 
acrocentric chromosome pair No. 12 in E. bruneus (white hollow arrows). Heterochromatin C-bands was 
consistent with the positive Ag-NORs sites on chromosome pair No. 2 in both fish species (thick black 
arrows). Heterochromatin C-bands were indiscernible for three pairs of chromosomes (Red arrows). 
Thin black arrows represent the chromosomes No. 24 in length. Scale bar = 5 [im. 
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Figure 5 Distribution of 1 8S rNDA on chromosomes of E. bruneus and E. moara by FISH. Distribution 
of 18S rNDA on chromosomes of E. bruneus (A) and E. moara (B) by FISH. White arrows indicate the 
biarmed chromosome, pairs No. 2, No. 9 and No. 24 in length, in both fish species. Scale bar = 5 ^im. 



percomorpha were mostly reconstructed the phylogenetic relationship by partial sequences 
of 18S rDNA but not 5S rDNA (data not shown). 

FISH analyses 

Multiple sites of 18S rDNA by FISH confirmed the data obtained by Ag-staining for 
NORs. In E. bruneus, six positive signals (both active and non-active) were identified, 
corresponding to metacentric (m) and submetacentric (sm) chromosomes (Fig. 5A). Four 
hybridization signals were observed on the short arms of submetacentric chromosomes in 
E. moara (Fig. 5B). 5S rDNA and 18S rDNA were found on different chromosomes. Two 5S 
rDNA sites were located on the arms of a medium-sized acrocentric chromosome pair in 
both£. bruneus and£. moara (Figs. 6A and 6B). 

Telomeric repeats of (TTAGGG) n showed the typically telomeric signals on both 
telomeres and/or centromeric region of all chromosomes in E. bruneus and E. moara. 
No positive signal was detected at interstitial sites (Figs. 7A and 7B). Ten chromosome pairs 
of E. bruneus were significantly stronger than the signals of the others (Fig. 7A). However, 
E. moara were characterized by uniform telomeric signals in strength and size (Fig. 7B). 

DISCUSSION 

Our previous study has distinguished E.bruneus and E. moara as two species based on 
their morphometric and skeletal characteristics {Guo et ah, 2008). The species-specific 
differences showed obviously on the bars of the body and stable skeleton characteristics of 
adult. And the pyloric caeca indicates their different feeding habits and digestive function, 
which means they could have different ecological niches. Mitogenome and molecular 
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Figure 6 Localization of 5S rNDA on chromosomes of E. bruneus and E. moara by FISH. Localization 
of 5S rNDA on chromosomes of E. bruneus (A) and E. moara (B) by FISH. White arrows indicate the two 
5S rDNA clusters located on the arms of one of acrocentric chromosome pair. Scale bar = 5 |im. 



comparisons confirmed E. moara to be a valid species of the family Epinephelidae 
(Liu et ah, 2013). Further, we developed a molecular method to differentiate both fish 
species (Guo et at, 2009). However, the cytogenetic backgrounds and evolutionary 
situation, which is very important to the cultivation and protection of fish resources, 
remains unclear for E. bruneus and E. moara. We here comprehensively analyzed the 
cytogenetic backgrounds, and reconstructed their phylogenetic relationships using 
molecular cytogenetic analyses. 

For E. bruneus and E. moara, species-specific characteristics presented in karyotype, 
NORs, C-banding and telomere distribution patterns. Karyotype variation appears to 
parallel speciation events in many groups of vertebrates (Morescalchi et at, 2007; Ruiz- 
Herrera, Farre & Robinson, 2012). Variations of NOR constituted a strong cytotaxonomic 
character in fishes (Fujiwara et al, 1998; Galetti, Aguilar & Molina, 2000). Many species in 
genus Epinephelus showed the same karyotypic characteristics, such as karyotype formula 
and NORs (Wang et al, 2012). However, karyotype formula and NORs were different 
between E. bruneus and E. moara (Table 1). Further, interstitial C-bandings were observed 
in E. bruneus, but not in E. moara. Similar interstitial heterochromatin was also found in 
E. coioides (Wang et al., 2010) and Diplectrum radiale (de Aguilar & Galetti, 1997). The 
distribution patterns of (TTAGGG)„ were different obviously between E. bruneus and 
E. moara. Cytogenetic differences were inter-specific, because E. bruneus and E. moara 
showed a similar geographical distribution (Heemstra & Randall, 1993; Guo etal, 2009), 
and coupled with heterogeneously morphological characteristics (Guo et al, 2008) and 
chromosomal structure. 
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Figure 7 Distribution of telomeric (TTAGGG)n sequence on chromosomes of E. bruneus and 
E . moara by FISH. Distribution of telomeric (TTAGGG)n sequence on chromosomes of E. bruneus (A) 
and E. moara (B) with telomeric (TTAGGG)n sequence using FISH. Red arrows indicate chromosomes 
with significantly stronger and larger telomeric signals than others in E. bruneus. Scale bar = 5 \im. 



The cytogenetic analyses suggested that the lineages of E. bruneus and E. moara recently 
derived within the genus Epinephelus. Both fish species share a uniform number of 
chromosomes to other species in the genus Epinephelus (Wang et at, 2012). However, 
they contained more biarmed chromosomes such as metacentric and/or submetacentric 
chromosomes (Table 1). In fishes, 48 uni-armed chromosome types like acrocentric 
chromosomes represented the ancestral complement of diploid origin (Ohno, 1974; 
Vitturi et at, 1991; Sola et at, 2000). In addition, most species in genus Epinephelus 
showed a conserved, NOR-bearing chromosome pair No. 24 (Table 1), while E. bruneus 
and E. moara showed additional NORs on chromosome pairs. For most vertebrates, the 
presence of a single NOR pair seems to be an ancestral character state (Hsu & Pardue, 1975; 
Schmid, 1978; Galetti, Molina & Affonso, 2006). Both E. bruneus and E. moara show even 
more constitutive heterochromatin (related to chromosomal rearrangements or variation) 
than other species in Epinephelus (Sola et at, 2000; Molina, Maia-Lima & Affonso, 2002; 
Phillips &Rab, 2001; Wang et at, 2012). 

Moreover, E. moara exhibited more plesiomorphic features than E. bruneus. The 
fundamental number (FN) of E. bruneus is larger than that of E. moara. Species with 
a larger FN are more derived in evolutionary terms (Martinez et at, 1989; Ghigliotti 
et at, 2007). Chromosomal rearrangements and genomic modifications were more 
obviously in E. bruneus compared to those in E. moara. Interstitial C-bandings appeared 
in E. bruneus imply the karyotypic rearrangement (Galetti, Aguilar & Molina, 2000), 



Guo et al. (2014), PeerJ, DOI 10.7717/peerj.412 



9/15 



PeerJ 



Table 1 Available cytogenetic data of the genus Epinephelus. 



Species 


In 


Karyotype formula 


FN 


NORs 


C-banding 


Reference 


E. adscencionis 


48 


48a 


48 


SCR(24) a 


C(l-24) 


{Molina, Maia-Lima & Affonso, 2002) 


E. akaara 


48 


5st + 43a 


48 


/ 


/ 


(Wanget al, 2004) 


ry „7 3 ■ 

b. alexandnnus 


48 


48a 


48 


bCR(24) 


NC(1-23J, aCR(24) 


(Martinez et al, 1989) 


E. awoara 


48 


48a 


48 


SCR(24) 


NC(l-23), SCR(24) a 


(Wanget al, 2012) 


■ ■ i 

b. bruneus 


48 


2m + 4sm + 42a 


54 


5>CR(24, 9, 2) 


NC (?) 
C(?) 

5>CR(24,9,2) 
SA(2) 
IR (?) 


Present study 


E. caninus 


48 


48a 


48 


SCR(24) 


/ 


{Rodriguez-daga, Amores & Thode, 

LyyJ ) 


C. LVIUlCltZb 


4ft 


1 cm _1_ 4Aa 








y vvu/ig ei ut., zuiu j 


E. diacanthus 


48 


2sm + 46a 


50 


i 


i 


(Natarajan & Subrahmanyan, 1974) 




48 


A m — 1- cm — 1- Act — 1- ^.d-n 


62 


i 


i 


yZ-idCflVy L*m KIT L*l, ZUUu ) 


E. fasciatomaculosus 


48 


48a 


48 


SCR(24) 


/ 


(Li & Peng, 1994) 


E. fasciatus 


48 


48a 


48 


SCR(24) 


/ 


(Li & Peng, 1994} 


H ti i cmcru ft nil ic 
IZ. j LlbLUVLlLLClLUj 


48 


_mii ^ud 


50 






( J inn pi nl 0OCif?\ 




4ft 


'iOd 


4ft 






ylvlllT llrltZz, CL ill., Lyoy) 


E. guttatus 


48 


48a 


48 


1 


/ 


(Medrano et al, 1988) 


E. lanceolatus 


48 


4st + 44a 


48 


1 


/ 


{Wanget al, 2003) 


F* 1 tr-i nln\,tiY,ri /c 
JZ. IflUlUCUl ILUi 


48 


48a 


48 


■ \p) 


C(\ -941 


\z^ua, in Zjiiuii, zuuj) 


E. marginutus 


48 


48a 


48 


SCR(24), TR(2) 


C(l-24) 

SCR(24) J 

TR(2) 


{Sola et al, 2000) 


E. merra 


48 


4m + 6sm + 4st + 34a 


62 


/ 


/ 


{Zheng, Liu & Li, 2005) 


E. moara 


48 


4sm + 44a 


52 


SCR(24, 9) 


NC(?) 

C(?) 

TR(?) 


Present study 


E. sexfasciatus 


48 


2sm + 46a 


50 


/ 


/ 


{Chen etal, 1990) 


E. tauvina 


48 


2sm + 46a 


50 


/ 


/ 


{Raghumath & Prasad, 1980) 



Notes. 

2n, diploid number; a, acrocentrics; C, centromeric; EA, nearly the entire arm; FN, fundamental number; IR, interstitial region; m, metacentrics; NC, almost 
indiscernible; NORs, nucleolar organizer regions; SA, short arm; sm, submetacentrics; st, subtelocentrics; SCR, subcentromeric region; TR, telomeric region; /, not 
available; ?, not mentioned or measured; Numbers in parentheses, the number of chromosome pairs. 
a Data estimated from illustrations and text in the respective papers. 



robertsonian rearrangements and/or reciprocal translocations (Eler et al, 2007). Despite 
the conservation of (TTAGGG)n sequence and location, slight changes in the telomeric 
sequences have occurred during vertebrate evolution {Meyne et al., 1990). Uniform 
telomeric distribution in E. moara is similar to other species in Epinephelus (Table 1 ) 
{Sola et al, 2000; Wang et al, 2012). However, remarkably high repetitions of telomere 
sequences seem to exist on ten chromosome pairs with stronger signals in 
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E. bruneus, which appear to involve complex homologous or/and non-homologous 
recombination. 

18S rDNA could be simultaneously applicable in the taxonomic and evolutionary 
analyses of groupers. The 5S rDNA seems to be unsuitable in the phylogenetic resolution, 
because the order percomorpha in the NJ trees were not recovered as monophylum. 
Distribution patterns of 5S rDNA of E. bruneus and E. moara are similar to other species, 
while that of 18S rDNA were different among species in genus Epinephelus (Sola et ah, 
2000; Wang et ah, 2010; Wang et ah, 2012). In addition, the different distributions of highly 
conserved 18S rDNA and telomere suggest the distinct genomes and evolutionary situation 
of the closely related species E. bruneus and E. moara. 

In summary, many useful cytogenetic charateristics are available to distinguish 
E. bruneus from E. moara, such as karyotypes, NORs, C-banding, 18S rDNA and telomere 
(TTAGGG) n distribution patterns. Moreover, the lineage of E. bruneus and£. moara seems 
to be derived recently, and E. moara exhibits more plesiomorphic features than E. bruneus. 
Molecular cytogenetic analyses could be applicable in identification of closely related 
species and reconstruct their phylogenetic relationships in groupers. 
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